Quantum metamaterial without local control 
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A quantum metamaterial can be implemented as a quantum coherent ID array of qubits placed in 
a transmission line. The properties of quantum metamaterials are determined by the local quantum 
state of the system. Here we show that a spatially-periodic quantum state of such a system can be 
realized without direct control of the constituent qubits, by their interaction with the initializing 
("priming") pulses sent through the system in opposite directions. The properties of the resulting 
quantum photonic crystal are determined by the choice of the priming pulses. This proposal can be 
readily generalized to other implementations of quantum metamaterials. 

PACS numbers: 81.05.Xj,74.81.Fa 
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Since the pioneering papers the investigation of 
metamaterials (quasicontinuous media built of artificial 
unit elements, that modify the properties of propagating 
electromagnetic waves) was driven by their unusual op- 
tical properties^ in a wide range of frequencies. They 
are used to produce antireflection coating and optical de- 
vices such as adaptive lenses, reconfigurable mirrors, con- 
verters, etc. Another direction of research is related to 
the development of metamaterials based on active units, 
e.g., Josephson junctions —ii. The emission of terahertz 
radiation due to the ac Josephson effect allows to build 
Josephson-based emitters, filters, detectors and waveg- 
uides operating in this spectral range important for appli- 
cations The nonlinear effective inductance of Joseph- 
son junctions makes them a convenient basis for nonlinear 
metamaterials — . 

Another fundamentally novel class of artificial media 
are quantum metamaterials, i.e., optical media, which 
maintain global quantum coherence over times exceed- 
ing the signal transition time and allow local control over 
quantum states and basic properties of their constituent 
elements^. The term reflects the fact that, similarly to 
classical metamaterials, these systems allow additional 
ways to control the propagation of electromagnetic fields, 
not available to standard materials. The optical proper- 
ties of the material are determined by its controllable 
coherent quantum dynamics. 

Recently various schemes of qubit-array control have 
been proposedi2r21. Quantum metamaterials were intro- 
duced ii>ii as an array of superconducting charge qubits 
placed inside a transmission line; the extensions to other 
implementations, including those in optical range, soon 
foUowedi^"— . In particular, experiments on a flux qubit 
inside a transmission line reproduce a number of atomic 
spectroscopy effects^^. Aside from the specific properties 
following from quantum metamaterial being an extended 
quantum system, it is also a medium with an optical re- 
sponse, which can be changed at will without changing its 
microstructure; such media are a focus of strong research 




FIG. 1: (Color online). ID quantum metamaterial in a su- 
perconducting transmission line. The superconducting is- 
lands (S) sandwiched between the superconducting strips 
form charge qubits (1 < n < A'^); D is the distance between 
the strips (S), Lo is the distance between the qubits. In the 
passive (in- and -out-) regions, with n < and n > N respec- 
tively, the transmission line parameters are chosen to optimize 
impedance matching with the active (metamaterial) section. 
The electromagnetic pulse propagates along z axis, A is the 
vector potential, H is the magnetic field of the electromag- 
netic wave. 



In this paper we show that the local manipulation of 
the quantum states of the constituent units of a quan- 
tum metamaterial is not the only possible way of ini- 
tializing it in a desirable state. To be specific, we show 
that a spatially-periodic quantum state of a ID quan- 
tum metamaterial (i.e., quantum photonic crystalii) can 
be reahzed by sending, from opposite directions, a pair 
of priming pulses through the system. This provides an 
easier way to initialize a quantum metamaterial, since 
it would not require the local control of each qubit in 
it. Disposing of this requirement significantly simplifies 
the experimental realization of quantum metamaterials 
by reducing the following: the complexity of their de- 
sign, the coupling to the environment, and therefore the 
decoherence due to both internal and external sources. 
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We model a ID quantum metamaterial using super- 
conducting charge qubits^^j, though the main idea of our 
approach is independent of the particular implementa- 
tion. The (identical) charge qubits are formed by super- 
conducting islands sandwiched between superconducting 
strips and separated from them by tunneling barriers (see 
Fig. HI). The metamaterial occupies the central ("active" ) 
section of the transmission line formed by the strips; in 
the rest of the system ("passive") qubits are replaced 
by capacitive elements chosen to optimize the impedance 



matching between the active and passive sections of the 
line. 

As in^^, we treat the electromagnetic field in the sys- 
tem classically, while quantizing the qubit degrees of free- 
dom. The field propagates along the z direction. Ne- 
glecting edge effects, the magnetic field will have only a 
y-component, and the vector potential in a section be- 
tween the nth and (n -t- l)st qubits can be chosen as 
■^xn&x- The energy of the active part of the system can 
be then written as 



Esc = 



Ej cos Lpn + 



+ cos Lp. 
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Here ipn is the superconducting phase of the nth island, 
Ej = $o^c/27rc and Uj = 2eIc/fiC are the Josephson 
energy and Josephson frequency respectively, Ic is the 
critical current, C is the Josephson junction capacity, 
$0 = hc/2e is the magnetic flux quantum, D is the dis- 
tance between the superconducting lines, W is the thick- 
ness of the superconducting strips. The first term in 
Eq. ([T]) is the charging energy of the nth qubit, the sec- 
ond is the Josephson energy, and the last two terms de- 
scribe the electromagnetic field energy. The time deriva- 
tive is indicated by a dot {A, cp) . We assume here that the 
qubit size is much smaller than the distance between the 
neighbouring qubits and that the vector potential does 
not depend on the coordinates in the area between them. 
Similarly, in the passive regions 
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Comparing the expressions for the electromagnetic field 
energy for the active [Eq. ^] and the passive [Eq. ([2|)] 
regions one can conclude that the impedance matching 
is fulfilled if C ie^Ej/fi^ijj. 

In order to obtain the Hamiltonian from the full clas- 
sical energy of the system. Esc + Epfc, we quantize the 
qubit contributions by writing 



Hr, 



{hujjf ( d 



Ej \dipn 



Hn |n) = En \n) , 

'2EjC0iHpn+hn. (3) 



Here the term hn describes the external controls of the 
state of the n-th qubit; its precise form is irrelevant here. 



Then we restrict the operator Hn to the subspace of its 
two lowest states, |0„) and |1„), while i?i^„ — £^o,ri — 
hsn- We can assume this due to the nonlinearity of the 
Josephson potential, assuming that both the temperature 
and the field amplitudes are small enough to keep only 
these states populated. 

As follows from Eq. ([1]), the interaction of the qubit 
with the electromagnetic field is given by 



Vn — 2Ej{\ — cosa„) COSipn 



(4) 



where a„ = 2TTDAxn/^o is the dimensionless vector po- 
tential. 

In the following, we consider the simplest case of fac- 
torizable wave function of the metamaterial. Each qubit 
is described by its own wave function \ipn)'- 

|,/'„(t)) =co(n,t)|0)e-*/2+ci(n,<)|l)e-*^*/2; (5) 

we omit here the subscripts in |0) and since the qubits 
are assumed to be identical. The Shrodinger equation 
for lipn) is thus reduced to a set of equations for the 
coefficients cq and ci: 

ihc^{n,t) = J2 (a|K(i) 1/3) c/^(n,i)e^("° (6) 

^=0,1 

where (a,/?) = 0,1, wi — e/2 and ujq = —e/2. From 
Eq. ([ij follows the equation for the electromagnetic field 
in the active region 1 < n < iV: 



a„ - v'^{an+i + a„_i - 2a„) + 
rsina„ (V'„|cos<^„ {ipn) = 0, 



where 



i2n^LaDEj ' 



32tt^c'^DE,, 



(7) 



(8) 
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In the passive regions we have instead 

a„ - u^(a„+i + a„_i - 2a„) = 0. (9) 

With the proper choice of parameters in the passive re- 
gion one can obtain u = v. 

For an analytical treatment, we make additional plau- 
sible assumptions. First, the electromagnetic field must 
be weak, a„ <^ I (otherwise the two-level approximation 
for qubits may not be valid); therefore 



Vn « Eja„ cos ipn 



(10) 



and sina„ « a„. Second, the field wavelength under a re- 
alistic choice of parameters greatly exceeds the interqubit 
distance and therefore allows to introduce a continuous 
variable z instead of the discrete label n (Ref. llj). Then 
Eq. ([7]) is transformed into: 



- v-^r^ + xiz,t)a = 0, 



(11) 



where v = vLq, a = a{z,t), and the "susceptibility" 

x{z,t)^r{^{z,t)\cosip\ij{z,t)) (12) 

is determined by the quantum state of the metamaterial. 

In agreement with our weak-field assumption, we an- 
alyze Eq. (HI]) in a quasilinear approximation. These 
analytical results will help illuminate the numerical solu- 
tions of Eqs. (|6|7I9|) for the following choice of parame- 
ters: Lq = 5 • 10^"' m, D ^ W ^ 10~^ TO, critical current 
/c = 4 • 10~^ A, e/2Tr - 5 • 10^° and Ej/hojj ^ 4, 
so, u ~ c. 

As shown in Ref. [ll| . the transmission properties of 
a ID metamaterial with a spatially-periodic quantum 
state, \^{z + Lm,t = 0)) = \i^{z,t = 0)), are those of 
a generalized photonic crystal. We therefore investigate 
the possibility of initializing the metamaterial without 
applying local controls to the constituent qubits. Specif- 
ically, let us initialize the metamaterial in the ground 
state (which can be done by, e.g., cooling) and send two 
priming pulses through the active region of the system in 
opposite directions. One would expect that due to the in- 
terference of these pulses and their action on the qubits, 
a spatially-periodic state of the quantum metamaterial 
would arise. For superconducting qubits the decoherence 
times are currently in the range 10-100 /is and the inter- 
level distance is ^ 10 GHz; the transition time for the 
pulse across a 1000-unit metamaterial (of length ~ 10^ 
cm) is ~ 3 • 10"^ s. Therefore, in the following, we can 
neglect decoherence effects. 

The backaction of the electromagnetic field on the 
qubits is quadratic in the field amplitude (Eg. llOp . there- 
fore the basic frequency of the priming pulses should sat- 
isfy the condition 2uj w e (see, e.g.^ii^li^). We assume 



that the spatial width of the pulse satisfies 
I > 2t:v/uj w a. 



(13) 



which ensures that each qubit in a superposition state 
undergoes many cycles of quantum beats while the pulse 
propagates past it. Here A is the wavelength of the prim- 
ing pulses. In this case, the metamaterial state evolution 
in the presence of the field can be investigated in the 
resonance approximation^^. So, Eqs. ([6|) become: 



ihco{z, t) = a^{z, t) (dooco + ^oicie '^*) , 
ihci{z, t) — a^{z, t) (dioCoe*^* -f duCi) , 



(14) 



where daf) — Ej {a\ cos(/?„ |/3) and z — hLq. In the case 
of sufficiently strong pulses the solution of Eg. dTT]) may 
be written as 



a^"'"' (z, t) = exp 



a*-^-' (z, t) = exp 



[z - {uj/k)tY 



{u/k)tf 



(15) 



where A is the amplitude of the pulses (the pulses have 
equal amplitudes), and (f>o is the initial phase. In order 
for our treatment of the metamaterial as a continuous 
medium to be consistent, the wave vectors should satisfy 
the condition A = 27r/fc ^ Lq, that is 



Lq <^ 2tiv/uj. 



(16) 



For V K, c and e/27r k, 5 • 10^° s^^, this yields Lq cm, 
which is a feasible requirement. 

Let all the qubits be initialized at i = in the ground 
state [co(z, 0) = l,ci(z,0) = 0], then, at later times, the 
quasi-monochromatic approximation [Eq. (jlSp ] gives 



\ci{z,t)\ = 



(^)|sin(/| 



n{z) 



■^{zflAt 



7(z) = A + 4A 



2 "00 



dii 



[cos(2fcz 



(17) 



(18) 



where A = 2w — e is the detuning from the resonance. 
The local Rabi frequency is 



\n{z)\ = 



h 



[cos(2A;z 



(19) 



The periodicity of 7(2) and Vl{z) [Eqs. [THl [TO] implies a 
spatially-periodic probability of exciting the qubit with 
period A/2. 

Figures [2] and [3] show the results of a numerical sim- 
ulation of such a process from Eqs. (|7l9p . Two prim- 
ing pulses at t = are located in the passive regions of 
the waveguide (Fig. [2]). After passing through the active 
region, they produce a periodically-modulated popula- 
tion of the levels of the qubits (Fig. |3]). A small distor- 
tion of the periodicity is due to the fact that the waves 
are non-monochromatic and can be improved by using 
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FIG. 2: (Color online). The priming or initializing pulses 
located in the passive regions of the waveguide line at the 
initial moment of time. The active region is in the center 
of the line between the two vertical dotted lines. The pulse 
parameters are v = c, uj = e/2, A = 0.18e, k = (l/25)-27r/Lo. 
The pulse width is / = 240I/o, and the pulse amplitude A = 
0.18. The matrix elements are doo = 0.4&, dn = 3.6&, and 
doi = Q.2he. The total length of the system simulated here is 
L = 2048Lo. 




FIG. 3: (Color online). Periodically- modulated population 
of the excited levels of qubits after passing the priming or 
initializing pulses through the system. The period of modu- 
lation, 131/0 ~ A/2, shows that this is due to the interference 
between these pulses. The slight aperiodicity is caused by the 
finite width of the pulse. 



wider pulses. For the system parameters chosen, the 
period of modulation in Fig. [3] is approximately equal 
to 13Loi in good agreement with the modulation period 
tt/A; = 12.5Lo following from the approximate solution 
Eq. The modulation amplitude can also be con- 

trolled by the pulse amplitude and the pulse width [cf. 
Eqs. (HZHHl)]. 

The periodic spatial modulation of the wave function 



\'(p{z)), created by the priming pulses, will affect the prop- 
agation of the subsequent probe pulse through the func- 
tion x(z) [Eq. ([T^ ] , which enters the wave equation pT|) : 



X{z,t) =xq + x[^ + cos(27rz/ L„J] + 



(20) 



where Lm is the modulation period, and xo depends on 
the quantum states of the qubits with minimal |ci|: 



Xo 



r 



■00 I Co, mill I 



rfii|ci, 



(21) 



For the case shown in Fig. [3l xo = {r/Ej)doo, since 
the qubits at positions such that cos(2fcz + (po) = — 1 
(i.e., where the local Rabi frequency is zero) remain in 
the ground state, and thus Ici^minl = 0. The quantity 
X describes the excitation of qubits at positions where 
cos(2fcz -I- 0o) = 1, and |ci| = |ci,inax|: 

X = r {dll - doo) (|ci,maxP - |ci,min|^) /2Ej. (22) 

For the case of Fig. [3] we have x = 
r (dll - doo) I Cl max 
If 



[(r/i?,7)doicS(z)ci(z)e^2Wi™]rf^| 



then one can neglect the third term in Eq. (f20|) to derive 
the dispersion equation and thus consider x as indepen- 
dent of time. So, x in Eq. (|20[) can be approximated 
as 



X{z) = Xo + X [1 + cos(27rz/L„)] 



(23) 



We seek the solution of Eg. pTjl with the above x(^) 
for small wave amplitude in the form of a Bloch wave^S: 



ak{z, t) = Uk{z) exp [i{kz - tot)] , 



(24) 



where Uk{z) is a periodic function with period Lm [i.e., 
Uk{z) = Uk{z -{- L„i)\. The dispersion equation is readily 
found: 



where 



Ln 



+ im/2 



x{z) exp (i2fcz) dz, k 



TTn 



(25) 



(26) 



It is clear from Eq. that a gap opens in the fre- 
quency spectrum for each k = im/ Lm (if the correspond- 
ing Fourier component of x(^) is different from zero, see 
FigH: 
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FIG. 4: Normalized dispersion relation aj(fc) for the quan- 
tum photonic crystal. Due to the space periodicity of x{^)t 
gaps open in the spectrum. Their positions and magnitude 
depend on -/-/m and X- This allows to control the parame- 
ters of the quantum photonic crystal by varying the quantum 
states of the qubits. The graph is plotted for Lm ~ 25Lo/2, 
doo = OAhe, dii = 3.6fc, doi = 0.2he, and v = c. The qubit 
population is the same as in Fig. [3] 



The gap appears for fc„=i — tt/ Lm, close to the frequency 

= '^{vkif + Wo. According to Eq. Wk„^^ = 

x/2, Wo = xo + X, and the gap equals 



Suji « x/,J{Trv/L„y + Wo. 



(27) 

In Fig. [S] the propagation of a probe pulse through a 
metamaterial with periodic xi^) is shown. One can see 
that while the pulse with the carrying frequency inside 
the bandgap undergoes significant reflection [Fig. [Slja)], 
this is not the case for a pulse with a frequency above 
the bandgap [Fig. (Sfb)]. 

In conclusion, we have shown that a ID quantum meta- 
material comprised of superconducting charge qubits in a 
transmission line can be initialized in a spatially-periodic 
state without excercising local control of the quantum 
state of individual qubits, by simultaneously passing 
through it two priming pulses in opposite directions. The 
modulation period is close to half the priming pulse wave- 
length, A/2. The subsequent probe pulse propagation 
through the resulting periodic structure was shown to be 
affected by the arising bandgaps. 

It is possible that such ID quantum metamaterials 
could be realized using the microwave transmission line 
that was recently experimentally investigated in^S, where 
a set of resonators coupled with qubits was connected 
with this line. If the transmission line is long enough, 
we can create a more complex picture of the qubit state 
distribution by varying the shape of the pump pulses. 
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